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1 Abstract 

2 It is well known that extreme or prolonged rainfall is the dominant trigger of landslides 

3 worldwide. While research has evaluated the spatiotemporal distribution of extreme rainfall and 

4 landslides at local or regional scales using in situ data, few studies have mapped rainfall- 

5 triggered landslide distribution globally due to the dearth of landslide data and consistent 

6 precipitation information. This study uses a newly developed Global Landslide Catalog (GLC) 

7 and a 13-year satellite-based precipitation record from TRMM data. For the first time, these two 

8 unique products provide the foundation to quantitatively evaulate the co-occurence of 

9 precipitation and landslides globally. Evaluation of the GLC indicates that 2010 had a large 

10 number of high-impact landslide events relative to previous years. This study considers how 

11 variations in extreme and prolonged satellite-based rainfall are related to the distribution of 

12 landslides over the same time scales for three active landslide areas: Central America, the 

13 Himalayan Arc, and central-eastern China. Several test statistics confirm that TRMM rainfall 

14 generally scales with the observed increase in landslide reports and fatal events for 2010 and 

15 previous years over each region. These findings suggest that the co-occurrence of satellite 

16 precipitation and landslide reports may serve as a valuable indicator for characterizing the 

17 spatiotemporal distribution of landslide-prone areas in order to establish a global rainfall- 

18 triggered landslide climatology. This study characterizes the variability of satellite precipitation 

19 data and reported landslide activity at the globaly scale in order to improve landslide cataloging, 

20 forecasting and quantify potential triggering sources at daily, monthly and yearly time scales. 
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1 1. Introduction 

2 

3 It is well established that intense or prolonged rainfall can trigger slope movements (Cannon and 

4 Ellen 1985; Caine 1980; Croizer 1986). These processes predominately occur within steep 

5 topography where intense or prolonged rainfall increases pore water pressures and decreases soil 

6 cohesion in the subsurface, causing the driving forces to overcome resisting forces on a hillslope 

7 and activate a landslide (Wieczorek 1996; Iverson 2000). Understanding the distribution of mass 

8 movement processes can be challenging as physically-based models require in situ knowledge of 

9 the surface and subsurface conditions at local scales in order to quantify how rainfall intensity 

10 and infiltration may trigger landslide events. In lieu of detailed surface information, research has 

11 relied on statistical or empirical comparisons of rainfall events and landslides to characterize the 

12 spatial and temporal distributions of mass movements at local or regional scales based on 

13 historical landslides and gauge-based rainfall (Caine 1980; Larsen and Simon 1993; Guzzetti et 

14 al. 2008; Lepore et al. 2011). A challenge inherent in both physical and empirical in situ 

15 evaluations is the availability of consistent precipitation information and landslide event data to 

16 effectively characterize the spatiotemporal distribution of landslide occurrences as well as 

17 validate these models, particularly over regional or global scales. 

18 

19 A newly developed Global Landslide Catalog (GLC) represents the first database of its kind to 

20 catalog ah rapidly-moving, rainfall-triggered landslides within the recent past at the global scale 

21 (Kirschbaum et al. 2009a). The catalog currently contains five complete years of data (2003, 

22 2007 - 2010) with continued reporting to the present. Through evaluation of this GLC dataset, 

23 we are able to extract information on the spatial and temporal frequency of landslide events at 

24 the global scale. While the GLC has several limitations identified below, the catalog provides a 
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1 foundation for exploring where and when landslide-triggering extreme storms have occurred 

2 over the globe and for characterizing hotspots for both landslides and extreme rainfall activity. 

3 We compare the GLC with satellite precipitation data from the Tropical Rainfall Measuring 

4 Mission (TRMM) Multi-Satellite Precipitation Analysis (TMPA), which offers a 13-year rainfall 

5 product at sub-daily time scales from 50°N to 50°S and spatial uniformity over most landslide- 

6 prone regions. 

7 

8 This study represents a first step in determining how variations in satellite precipitation are 

9 related to variations in landslides reported in the GLC, primarily on seasonal and inter-annual 

10 time scales and on regional spatial scales. Although satellite -based rainfall information has 

11 limitations in mountainous regions, the globally uniform nature of the data makes it very useful 

12 in comparison with the global landslide database. Previous research has applied remotely sensed 

13 data to evaluate magnitude-size relationships, progression and cataloging of landslides as well as 

14 to study channel morphology and contributing slope area to potential failure sources (Lashermes 

15 et al. 2007; Galewsky et al. 2006; Singhroy et al. 2002; Petley et al. 2002; Haeberlin et al. 2002). 

16 This analysis as well as previous related studies employ satellite-based rainfall estimates to 

17 evaluate landslide hazards with the goal of assessing their distribution over the global or regional 

18 scale (Hong et al. 2006; Kirschbaum et al. 2011; Liao et al. 201 1). While the spatial resolution 

19 (0.25°x0.25°) of the TMPA product precludes its use for detailed hillslope investigations since 

20 precipitation can vary substantially within a single grid box, TMPA data shows promise in 

21 characterizing landslide processes over larger areas using statistical or empirical methodologies. 

22 
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Upon completion of the rainfall-triggered landslide catalog for 2010 1 , the authors noted 
significantly more high-impact landslides for 2010 compared to previous years in the record. For 
example, a catastrophic mudslide occurred in Zhouqu County in Gansu, China on August 8 th , 
2010 which killed 1,765 people and resulted in an estimated 759 million USD in damages (EM- 
DAT 2011). Additional damaging landslides occurred in: Bududa, Uganda in March causing 
nearly 400 fatalities; Leh in Ladakh, Indian Kashmir in August which resulted in an estimated 
245 fatalities; and a series of events in eastern Brazil during January and April which killed over 
700 people. Media reports identified intense rainfall as the trigger for each of these events, which 
mobilized large volumes of material and interacted with the local morphology to generate 
catastrophic landslides. 

Drawing upon the global nature of the GLC as well as quasi-global satellite precipitation 
information, this research seeks to determine whether 2010 was an anomalous year for extreme 
precipitation and landslide activity as well as outline potential sources for this behavior. This 
research examines the co-occurrence of the GLC and TMPA precipitation in order to 
quantitatively determine how these datasets may inform each other in terms of the spatial 
distribution of extreme rainfall and occurrence of landslide “hotspots” over the globe. Through 
this evaluation, the analysis also considers how these established relationships may help to 
potentially forecast landslide activity and variability at seasonal, annual and decadal scales. This 
work may also serve as a building block to move one step closer to developing a global 
climatology of rainfall-triggered landslides, which currently does not exist and is greatly needed 
by many different organizations from international aid agencies to local governments. 

1 Landslide inventory information and documentation is available at 
http://trmm.gsfc.nasa.gov/publications dir/potential landslide.html 
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This paper focuses on observations of anomalous rainfall-triggered landslide reports during 2010 
and considers how the landslides relate to mean monthly or daily rainfall for 2010 over three 
particularly active areas. The paper then considers the extent to which corresponding extreme 
daily or monthly rainfall signatures differ in 2010 compared to previous years in the GLC. 

Lastly, this study provides a discussion of the potential sources for why 2010 may represent an 
anomalous year over the three study areas considered as well as how this type of analysis may be 
expanded and applied in the future. 

1.1. Data description 
Landslide inventory 

Few databases have attempted to catalog landslide occurrences at the global scale outside of 
merely listing the sources for relevant landslide articles. Petley et al. (2005) has developed a 
valuable global database of fatal landslide events from 2003 to the present and reports on recent 
significant landslides around the world on a blog site (http ://blo gs . agu . or g/landslideblo g/) . The 
Global Landslide Catalog (GLC), developed by the authors, considers all rapidly-moving 
landslides (term used herein to refer to debris flows, mudslides, landslides, etc.) directly 
triggered by intense or prolonged rainfall (Kirschbaum et al. 2009a). Landslide event information 
is obtained from online media reports, disaster databases, and governmental and non- 
governmental organizations, as well as personal correspondence in some cases. The landslide 
entries include information on the date of the landslide, the location (both nominal and 
latitude/longitude), type of movement (if available), trigger (heavy rainfall, storm name, or any 
secondary triggers if reported) and impacts (fatalities, injuries or affected persons, and additional 
information). The GLC has been compiled since 2007 and provides a retrospective assessment 
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1 for 2003 (Kirschbaum et al. 2009a). The GLC has also been used to evaluate a real-time, quasi- 

2 global estimation of landslide events using satellite precipitation information (Kirschbaum et al. 

3 2009b). 

4 

5 This inventory provides the first global picture of all available rainfall-triggered landslide 

6 reports; however, the catalog only represents a fraction of the total number of rainfall-triggered 

7 landslides occurring around the world due to several limitations. The primary challenges of this 

8 cataloging effort stem from the complex nature of landslide processes as well as the availability 

9 and accuracy of landslide reports. The catalog only includes a landslide report if rainfall was 

10 identified as the primary trigger of the event. The GLC relies on media reports and is 

11 consequently impacted by reporting issues including accuracy of the reported information and 

12 challenges in identifying the timing and location of reported events. While the inventory contains 

13 some information gleaned from non-English articles, the GLC primarily uses landslide reports in 

14 English. Landslide information and impacts are also frequently grouped with other hazards (e.g. 

15 floods, tropical cyclones), making it difficult to clearly identify the timing, location and 

16 magnitude of the specific landslide events. Lastly, it is often difficult to identify the precise 

17 location of the landslide event due to vague reporting or difficulty in locating remote villages 

18 where events have taken place. A qualitative “confidence radius” metric is included to indicate 

19 the relative confidence of each report’s latitude and longitude. An additional qualitative metric is 

20 used to describe the relative size of reported landslides based on reported impacts and the areal 

21 extent affected (i.e. street, town, or larger) with the goal of discriminating between smaller and 

22 larger events. 

23 
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1 Despite the cited challenges, the landslide catalog contains over 2,700 events with 10,500 

2 reported fatalities for 60 different countries over the years 2007 to 2010. Fig. 1 displays the 

3 number of landslides that caused at least 1 fatality (referred to herein as fatal landslides) per 

4 country over the consecutive GLC record as well as fatal and total landslide reports by month for 

5 each year. Upon compilation and evaluation of the 2010 record, we observed a notable increase 

6 in the number of reported landslides, fatal events and fatalities, including a three-fold increase in 

7 the number of reported events and a two-fold increase in reported fatalities and fatal events as 

8 compared to previous years. In this study we consider how reported and fatal landslide events for 

9 2010 and previous years may co-vary with extreme or prolonged rainfall in order to establish a 

10 potential indicator for more effectively characterizing landslide -prone areas at the global scale at 

11 seasonal and interannual time scales. 

12 

13 Rainfall Information 

14 This research uses daily TMPA precipitation data to characterize rainfall signatures and 

15 variability that produce damaging landslides. This merged satellite-based precipitation product 

16 provides a 13-year, 3-hourly continuous record from 50°N to 50°S at 0.25° x 0.25° resolution 

17 every 3-hours. The TMPA (Version 6) rainfall analysis uses multiple satellite estimates, all 

18 calibrated or adjusted by the TRMM radar/radiometer combined estimate (TRMM product 

19 number 2B31) and also uses a monthly rain gauge analysis to adjust the bias over land areas 

20 (Huffman et al. 2007, 2010). The TMPA has been validated against daily gauges and does well 

21 at reproducing the high end of the daily rainfall distribution. It also has the advantage of 

22 uniformity over the globe. The TMPA has also been shown to be useful in flood detection when 

23 it is used to drive a hydrological model (Yilmaz et al. 2010). Although there are limitations 
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1 regarding the accuracy of satellite rainfall estimates, including merged products such as the 

2 TMPA, the daily and monthly satellite -based estimates used here should be adequate for 

3 comparison with the landslide information. 

4 

5 The GLC was initially developed for evaluation of a global landslide hazard forecasting 

6 algorithm, which couples a global static landslide susceptibility map with TMPA satellite-based 

7 rainfall intensity and duration information to identify potential areas of landslide activity (Hong 

8 et al. 2007; Kirschbaum et al. 2009b; Hong et al. 2006). Hong et al. (2006) calculated an 

9 empirical intensity-duration (I-D) threshold using TMPA data and a set of global landslide 

10 events to specify an average rainfall intensity threshold above which a landslide may be 

11 triggered. This study uses the global 1-day threshold value of 79 mm/day to represent potential 

12 landslide triggering due to extreme rainfall. 

13 2. Rainfall Anomalies and 2010 Landslides 

14 2.1 Global Distribution of Landslides 

15 

16 In order to characterize the relationship between extreme precipitation hotspots and landslides 

17 during 2010, we first consider whether the pronounced increase in landslide reports (either fatal 

18 or total reported events) represents an artifact of the catalog or if there are observable patterns in 

19 increased activity for 2010. Fig. 1 highlights the increase in the number of reports, fatalities and 

20 fatal events for 2010. Fig. 2 displays the distribution of reported landslides and fatal landslides, 

21 plotted by month for the years 2007-2010 for the two large regional areas (South Asia and the 
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1 Americas) that dominate the statistics. The two figures illustrate that many of the landslide 

2 reports are distributed in reasonably well-defined regional clusters. 

3 

4 We identify three key areas where reported landslide activity has been fairly consistent 

5 throughout the record but which exhibit a pronounced increase in reports during 2010. These 

6 regions include Central America, the Himalayan arc and central-eastern China, representing 

7 some of the most active rainfall-triggered landslide areas in the world (Fig. 2). Through this 

8 evaluation, we investigate the connection between increased reporting and anomalous rainfall 

9 activity in these regions and how that may affect the global total of landslides. These areas are 

10 chosen based on the availability of landslide inventory information; however, we feel that the 

11 test areas provide a representative cross-section of highly susceptible areas over the globe and 

12 cover diverse climatologic and topographic regimes. Landslide reports and TMPA pixels were 

13 extracted for each of these regions and compared for both extreme daily rainfall and monthly 

14 anomalies. Several other regions displaying regional maxima and minima in 2010 were 

15 evaluated but are not included in this paper because they contained a limited number of data 

16 points for other years in the record or the landslide reporting was deemed to be inconsistent. 

17 Areas evaluated include the Northwest and Appalachian range within the United States, parts of 

18 South America, the Philippines, Indonesia and Southeast Asia. 

19 

20 2.2 Satellite Rainfall 

21 Fig. 3 illustrates the TMPA yearly rainfall anomaly map for 2010, using the yearly average over 

22 the evaluation period 1998-2010. Large anomalies over Burma and central Africa are primarily 

23 due to poor gauge coverage in the gauge analysis used by the TMPA in the current version (to be 
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1 corrected in the next version of this product). The climatology and anomalies of the TMPA 

2 compare favorably with exclusively gauge-based global products over the three test areas in this 

3 study. While the TMPA product has known problems over orographically complex terrain due to 

4 challenges in passive microwave rainfall retrievals, using the TMPA data allows for consistency 

5 when computing monthly totals and daily extreme precipitation statistics over different regions. 

6 

7 The global anomalies for 2010 show above-average rainfall over the three study areas (Central 

8 America, China, and the Himalayan arc) for 2010. However, to explore how the various time 

9 scales and rainfall intensities are related to landslides events, we examine the distribution of both 

10 monthly and extreme daily rainfall using three test metrics: 

11 

12 Monthly Rainfall Anomalies 

13 The 2010 and other years’ monthly rainfall totals were compared to the month climatology 

14 calculated from the TMPA 3-hr resolution record for 1998 - 2010 to obtain anomaly fields for 

15 comparison with the landslide data for 2010 and other years. 

16 

17 Daily threshold exceedance 

18 The 1-day rainfall intensity value from the global I-D threshold is used as the threshold to 

19 determine how frequently extreme rainfall occurred over the test areas. Any time the daily 

20 precipitation for a given pixel exceeds 79 mm/day, it is considered a “hit.” The number of “hits” 

21 are summed over the test area by month and divided by the number of total pixels in the test area 

22 to provide a relative threshold exceedance rate, which is intended to provide a comparison 

23 between extreme daily rainfall for 2010 and previous years. Exceedance rates are computed 
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1 monthly for 2010 and averaged for the years 2007 - 2009 to be consistent with the continuous 

2 GLC record. These values are compared with reported landslides over the same month for each 

3 region. The number of “hits” is also summed over each study region for each month and 

4 compared to monthly precipitation and fatal landslides. 

5 

6 Quantile -Quantile Plots 

7 The third metric tests whether 2010 daily rainfall values are statistically significantly different 

8 from previous years for the upper tail of distribution for the TMPA record. Precipitation 

9 quantiles are calculated for daily rainfall for 2010 and the years 1998 - 2009 within each study 

10 area and plotted on Quantile-Quantile (Q-Q) plots to determine if the probability distributions of 

11 two samples are independent. The two time periods are considered to be from different 

12 distributions if the quantile values diverge from their joint linear distribution. Quantiles are 

13 plotted against two lines: the 1:1 line (green) has a slope of 1, and interquartile line (red) shows 

14 the linear distribution of the 25 th and 75 th quantile for both datasets (shown in Figs. 4 - 6). A 

15 steeper positive slope of the interquartile line indicates that 25 th and 75 th quantiles of the 2010 

16 precipitation data have a larger spread (i.e. more extreme values). If the interquartile line 

17 diverges significantly from the 1:1 line, it suggests that the distributions between the two datasets 

18 (2010 vs. 12-year record) are different within the interquartile range of each dataset. 

19 

20 The Kolmogorov-Smirnov (K-S) test is then used to compare the probability distribution of the 

21 two datasets by calculating the distance between the cumulative distribution functions of the two 

22 samples (Massey 1951). The null hypothesis for the K-S test assumes that the two datasets come 

23 from the same continuous distribution. The null is rejected if the two datasets have different 
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1 continuous distributions at a given significance level (a) based on the K-S test statistic and p- 

2 value. The K-S test statistic is defined as the maximum difference between the two datasets’ 

3 cumulative distributions and the corresponding p-value determines the probability of obtaining 

4 the given K-S test statistic. The Q-Q distribution plotting and K-S test are performed using 

5 Matlab software. If the p-value is lower than the designated significance level, the null 

6 hypothesis is rejected. Since this research is focused on comparing only extreme daily 

7 precipitation, the K-S test is computed for precipitation values exceeding the 75 th quantile of the 

8 precipitation record. Q-Q plots and K-S test statistics are calculated for each region and shown in 

9 Table 1. 

10 3. Landslides and precipitation hotspots 

11 3.1 Central America 

12 The Central American test area extends from the southern tip of Mexico to Costa Rica and 

13 includes 355 TRMM (approximately 221,900 km 2 ) pixels and 86 landslides from 2007 - 2010 

14 (Fig. 2a, b). The monthly climatology shows a peak in boreal summer rainfall, punctuated by a 

15 mid-summer drought (MSD) in July, consistent with previous research (Magana et al. 1999) 

16 (Fig. 4a). Tropical cyclone activity is somewhat suppressed during the MSD and picks up again 

17 in late August or September. 

18 

19 The 79 mm/day minimum threshold was applied for the years 2007 - 2010 to evaluate daily 

20 exceedance values; however, the global threshold proved to be too high for the daily 

21 precipitation values observed in this area, resulting in only a few days when the threshold was 

22 exceeded. Recent work has suggested that a regionally-based I-D threshold may be better 

23 equipped to identify potential landslide triggering conditions over this study area, citing a value 


13 



1 of 39 mm/day as a more appropriate minimum daily rainfall threshold (Guzzetti et al. 2008; 

2 Kirschbaum et al. 201 1). Fig. 4b plots the rainfall threshold exceedance rate for 2010 and 2007 - 

3 2009 using the regional threshold proposed by Guzzetti et al. (2008) along with corresponding 

4 reported landslides. The 2010 exceedance rate highlights a dual peak in extreme precipitation 

5 that nearly parallels the occurrence of landslides reported in 2010. Both the exceedance rate and 

6 reported landslide values are nearly twice as large for most of the summer months in 2010 

7 compared to the same months in 2007 - 2009. 

8 

9 Fig. 4c plots the quantile values for the 12-year TMPA daily record (x-axis) and the 2010 daily 

10 precipitation values (y-axis). Table 1 provides results from the K-S test for values above 75 th 

11 quantile, showing a K-S test statistic of 0.1792 and p-value of 0.0026. These values suggest that 

12 the null hypothesis can be rejected at the 99.7% confidence level. From the results in Table 1 and 

13 Q-Q plot in Fig. 4c it appears that 2010 and the previous record have different distributions 

14 above the 75 th quantile and that the precipitation quantiles are somewhat larger for the 12-year 

15 record compared to 2010. However, because there is a significant positive difference (on the 

16 order of 3 mm or larger) between the interquartile range line (red) and the 1:1 line (green), the 

17 results indicate that the 2010 daily precipitation quantiles are actually larger when compared to 

18 the 1998 - 2009 record. 

19 

20 3.2 Himalayan Arc 

21 Along the southern margin of the Himalayan mountain range, including portions of India, Nepal, 

22 and Pakistan, monsoon rains trigger large numbers of damaging and fatal landslides each year. In 

23 2010, landslides caused approximately 500 fatalities in July through September over the study 
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1 region. The study area for this evaluation contains 700 TMPA pixels and 284 landslide reports, 

2 which covers an area of roughly 468,000 km“ (Fig. 2c, d). Fig. 5a displays the monthly 

3 climatology for this region and shows a clear 100-150 mm higher peak in monthly precipitation 

4 during July through September for 2010 compared to the climatology. Fig. 5b plots the daily 

5 threshold exceedance rates for 2010 and 2007 - 2009 using the global 79 mm/day threshold. 

6 When exceedance rate values are compared with the reported landslides, results show that 

7 exceedance rate values were approximately 1.5 times higher than the average values from 

8 previous years for the months of July through September. The number of reported landslides 

9 shows a similar peak, with values nearly five times higher for 2010 compared to the mean of 

10 previous years, and roughly twice as high for fatal landslides over the same time period. 

11 

12 The Q-Q plot shown in Fig. 5c indicates that quantile values for the 2010 data diverge from the 

13 interquartile line as well as the 1:1 line after approximately 3.7 mm/day, corresponding to the 

14 79 th quantile. Results from Table 1 indicate that the K-S test produces a high K-S test statistic 

15 and a very low p-value, suggesting that the null hypothesis may be rejected at the 99.9% 

16 confidence level and that the extreme precipitation values for 2010 are significantly higher than 

17 for the 1998 - 2009 TMPA record. 

18 

19 3.3 China 

20 The test area within central-eastern China contains 810 TRMM pixels (approximately 512,700 

21 km“) and 34 landslides (Fig. 2c, d). Fig. 6a displays a pronounced peak in the 2010 monthly totals 

22 for July and August, which is consistent with the peak in landslides during the same months (Fig. 

23 6b). The rainfall threshold exceedance rates for 2010 and 2007 - 2009 indicate that July is the 
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1 peak month for extreme daily precipitation. However, when comparing the monthly values with 

2 the landslide record, it is evident that anomalously high rainfall accumulations were observed for 

3 both July and August. The Q-Q plot shows that after approximately 8.7 mm/day (corresponding 

4 to the 95 th quantile) the 2010 quantile values diverge from the 12-year distribution, suggesting 

5 that the most intense daily precipitation values were higher in 2010 compared to previous years 

6 (Fig. 6c). The K-S test for the 75 th quantile and higher (corresponding to a rain rate of 2.42 

7 mm/day) does not reject the null hypothesis. However, at the highest precipitation values (above 

8 the 90 th quantile) the null is rejected with a p-value of 0.0276 at the 96% confidence level. While 

9 the climatology and highest daily precipitation values indicate that 2010 may be different from 

10 previous years, this area provides much less conclusive results. Sources of uncertainty are 

11 discussed below. 

12 3.4 Comparison of the three test regions 

13 Fig. 7 compares the monthly rainfall and exceedance threshold values for each month in the 

14 record over the three study areas. The 79 mm/day threshold was used for the India and China 

15 study regions and the 39 mm/day threshold was applied for Central America. Fig. 7a displays a 

16 scatter plot of monthly rainfall vs. exceedance values for each month over the study regions from 

17 2007 - 2010, showing a clear positive linear trend between increasing monthly rainfall totals and 

18 increased number of ‘hits’ when the daily rainfall threshold was exceeded. Monthly rainfall 

19 (Fig. 7b) and exceedance threshold values (Fig. 7c) are also compared to fatal landslides for each 

20 month. The number of fatal landslides is averaged over each 50 mm or 50 exceedance value 

21 interval. Fig. 7b and c suggest that despite having an uneven number of data points within each 

22 bin, there appears to be a slight increase in the average number of fatal landslides as the monthly 

23 rainfall or exceedance values increase. 
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1 4. Discussion 

2 The GLC dataset provides a unique global validation proxy for evaluating co-occurrence of 

3 extreme and prolonged rainfall and high-impact landslide events. Within this evaluation, we 

4 identify 2010 as an active year for rainfall-triggered landslides at the global scale and relate 

5 precipitation signatures to the GLC in order to determine how inter-annual precipitation 

6 variations are related to variations in landslides within the three identified landslide-prone 

7 regions. While it is well-known that intense or prolonged precipitation and landslide initiation 

8 processes are linked, the global nature of the GLC and TMPA precipitation record allows us to 

9 quantitatively diagnose this relationship at regional and global scales for the first time. 

10 Establishing direct correlations between these two products in terms of how they co-vary over 

11 space and time is complicated due to incomplete data records. This work will continue to 

12 improve upon the existing GLC in order to amass a more robust record of landslide events and 

13 accurately link rainfall patterns with landslide triggering events. 

14 

15 From this analysis, we determine that there is a clearly observable increase in rainfall-triggered 

16 landslide reports during 2010, compared to previous years. Fig. lb and c displays the monthly 

17 distribution of fatal and total reported landslide for the years 2007 to 2010. The increased peak in 

18 fatal reports during August corresponds to a large peak in activity from monsoon rains, with 

19 approximately 60% of the fatal reports occurring in China, Nepal and India alone. The total 

20 number of landslide reports is on the order of three times larger than the previous years’ 

21 inventory. In addition to the increase in anomalous rainfall-triggered events observed over the 

22 three study regions in 2010, the increased number of events may also partially be a result of 

23 improved reporting and better cataloging of reports. One way to consider a more realistic global 
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1 distribution of landslide activity is to only consider fatal reports, a potentially more reliable 

2 statistic since fatal events are generally more likely to be reported. Despite the short record, we 

3 have observed an increase in the number of fatal landslides over the lifetime of the inventory, 

4 which is consistent with Dr. Petley’s findings for fatal landslides for 2003 to 2010 and shows a 

5 peak in fatal landslides for 2010 (Petley 201 1). While variability in reporting accuracy is 

6 extremely challenging to characterize between regions, we anticipate that as we have more years 

7 of landslide report data we may be able to classify geographic biases in the GLC. 

8 

9 The precipitation anomalies shown in Fig. 3 highlight several areas that have experienced 

10 particularly wet seasons in 2010, including the three study areas evaluated here as well as 

11 southern India, Indonesia, eastern Australia, and northwestern South America. Within these 

12 areas, Indonesia, and portions of Colombia, Venezuela and Brazil in South America also 

13 experienced increases in landslide activity with more fatal landslide reports. Comparatively, in 

14 countries with negative precipitation anomalies such as Vietnam, we observed fewer reported 

15 landslide events. 

16 

17 There are many driving factors influencing regional variations in rainfall accumulation and 

18 intensity on seasonal and annual scales. The El Nino/Southern Oscillation (ENSO), while global 

19 in nature, has highly variable impacts on precipitation accumulation at regional scales 

20 (Ropelewski and Halpert 1987; Curtis and Adler 2003). The Nino 3, 3.4 and 4 indices show a 

21 large positive anomaly in January and February 2010, suggesting a strong El Nino (NOAA 

22 201 1). El Nino conditions continued until late February when the ENSO indices indicated a 

23 transition into La Nina conditions beginning in July and peaking in the mid to late fall, 2010. 
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1 Within the United States, wet weather likely amplified by El Nino and La Nina conditions and 

2 contributed to an increased number of landslide reports in southern California in January and 

3 February and California and Washington in December. While the El Nino signal was strongest at 

4 the beginning and the La Nina signal was strongest at the end of 2010, the majority of the rainfall 

5 events associated with anomalously high landslide activity over the study regions occurred in the 

6 boreal summer months, coinciding with a fairly weak ENSO signal. Below we discuss the impact 

7 of ENSO signals within the three study regions and their possible delayed impacts on boreal 

8 summer precipitation. 

9 

10 Central America 

11 The 2010 values during the summer months show a 50 - 100 mm increase in accumulation 

12 compared to the 12-year climatology, with the largest peak in August and September over this 

13 region. The landslide reports show a similar peak in reporting during May, August - September, 

14 and November for 2010 (Fig. 4b), with three times more fatal landslides and over four times 

15 more total reports. The extreme daily rainfall quantiles and monthly accumulations all suggest 

16 that the increase in reports tends to mirror the observed anomalous precipitation in the TMPA 

17 record. The peak in reports during May and November were likely the result of two tropical 

18 cyclones: tropical cyclone Agatha on May 29-30 th caused approximately 9 fatal landslides in 

19 Guatemala, and tropical storm Tomas in early November caused 2 fatal landslides in Costa Rica 

20 and many other landslide reports along roads. Due to the extreme nature of these events, there 

21 may also have been an over-reporting bias for these storms. 

22 
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1 One of the reasons for the positive precipitation anomalies over this region could result from a 

2 fairly active 2010 tropical cyclone season in the tropical Atlantic. ENSO has been shown to 

3 modulate interannual tropical cyclone frequency and redistribute precipitation extreme s(Elsner et 

4 al. 1999; Curtis et al. 2007). Curtis (2002) found that during the summer before a La Nina event, 

5 such was the case for 2010, precipitation follows a similar pattern to El Nino or Neutral patterns 

6 at the beginning of the summer, but then increases considerably in September. This inter-annual 

7 pattern can be linked to sea surface temperature changes and moisture due to ENSO as well as 

8 enhanced tropical cyclone activity. 

9 Himalayan Arc 

10 Monthly and extreme precipitation signals for the Himalayan Arc study area point to increased 

11 precipitation totals during the summer monsoon months in 2010, with the null hypothesis being 

12 rejected at the 99.9% confidence level and the exceedance values indicating a nearly two-fold 

13 increase in the number of extreme precipitation days during 2010 compared to the 2007-2009 

14 period (Fig. 5). Several studies have evaluated the connection between monsoon rains and 

15 landslide susceptibility over this region (Nagarajan et al. 2000; Gabet et al. 2004; Petley et al. 

16 2007). Indian Monsoon rainfall has been shown to strongly correlate with ENSO phases due to 

17 the coupling of tropical ocean-atmospheric modes over the Indian Ocean (Krishnamurthy and 

18 Goswami 2000). While Indian monsoon conditions are often suppressed during an El Nino event 

19 (Krishna Kumar et al. 2006; Webster et al. 1998), in the summer following a strong El Nino 

20 event, there is a tendency for above-normal precipitation with the most pronounced signal in 

21 August and September (Park et al. 2010). While ENSO is not the only circulation pattern 

22 contributing to the variability of boreal summer rainfall, results indicate that the strong ENSO 
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1 signal during 2010 may have played a sizeable role in the positive precipitation anomalies 

2 observed over this region. 

3 

4 China 

5 Results from Fig. 6 show that there is a pronounced peak in monthly rainfall during July and 

6 August, corresponding to an increased number of landslide reports. However, both exceedance 

7 values and rainfall quantiles do not clearly show the relationship between precipitation extremes 

8 and landslide reporting. The inconsistency in the landslide reporting as well as the size of the 

9 study area may be the limiting factors in this evaluation since only 34 landslides were reported 

10 over a very large area (5 12,700 km“) during the four year period. This is likely an 

11 underestimation of the GLC, due to reporting or language barriers, the occurrence of landslides 

12 in remote areas, or the influence of other triggers (such as the Wenchuan earthquake in 2008). In 

13 addition, triggers such as antecedent moisture or short, intense rainfall events (less than 24 

14 hours), such as what caused the Zhouqu mudslide, were either not included in the database or not 

15 adequately resolved by the satellite information. 

16 

17 Evaluating the sources of seasonal and annual variability of summer precipitation over China is 

18 challenging due to diverse climate zones and multiple ocean-atmosphere feedbacks influencing 

19 precipitation in this region. The Asian monsoon has been shown to strongly couple with tropical 

20 SSTs and the propagation of atmospheric circulation over the western Pacific, which affects the 

21 modulation of the Asian monsoon (Yang and Lau 2004). Directly north of the China study area, 

22 Feng and Hu (2004) one study also found that during a strong ENSO, there is a coupled 

23 relationship between the Indian summer monsoon and precipitation variability over northern 
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1 China. The strength of the precipitation-landslide signal in this region is complicated by a 

2 number of factors, including a limited number of landslide reports available over this region with 

3 varied reporting accuracy, the large size of the study area considered, and additional triggers 

4 influencing landslide susceptibility and triggering including anthropogenic impacts such as 

5 population density, road construction, and mining. 

6 

7 Each of the study areas is impacted by a different set of regional atmospheric circulation 

8 patterns, annual rainfall totals and surface susceptibility characteristics. Figs. 2, 4 - 7 and Table 1 

9 suggest that within the three regions, there is a statistically significant positive anomaly in 2010 

10 when comparing landslide reports and rainfall signatures within the TMPA record. It is clear 

11 from the analysis of these test areas that the total increase of 2010 landslides globally is related 

12 to changes in precipitation over different areas. However, the causes of these precipitation 

13 anomalies associated with increases in landslides vary from region to region. There is no 

14 singular reason for the precipitation increases and the related landslide maximum in 2010. 

15 

16 This analysis is a first step and should be considered preliminary for a number of reasons. First, 

17 because this evaluation only considers four years of data, identify temporal signals in landslide 

18 reporting may produce erroneous results over some regions. The increase in landslide reports 

19 over the period of record is observed by another database (Petley 201 1), suggesting that the a 

20 signal exists despite regional heterogeneities. As this rainfall-triggered landslide inventory 

21 continues to increase, it will provide more information to better quantify the regional reporting 

22 biases inherent in this type of a catalog. 

23 
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1 Second, merged satellite products offer a unique perspective on rainfall distribution by providing 

2 an inter-comparison framework amongst regions and through time. However, the sampling 

3 frequency of current microwave sensors does not allow for continuous monitoring of 

4 precipitation features and as a result, short events or peak intensities may not be accurately 

5 resolved by space-bome instruments or merged data sets. Comparing TMPA with the gauge- 

6 based products indicates that both products adequately highlight regional precipitation 

7 anomalies, but may not always resolve the exact magnitude of precipitation intensity. Comparing 

8 the relative magnitude of cumulative or daily exceedance values amongst regions allows for 

9 more consistent evaluation of the global prototype landslide algorithm system and evaluation of 

10 the rainfall-landslide relationship. This underscores the motivation for identifying an observable 

11 connection between the GLC and TMPA data so as to develop a potential indicator for high- 

12 intensity rainfall, particularly over mountainous regions where existing products may have 

13 difficulty accurately resolving precipitation. 

14 

15 Third, antecedent moisture may also play a sizeable role in the initiation and distribution of 

16 landslide events. Moisture within the soil can cause a buildup of pore water pressure such that 

17 smaller rainfall events occurring when the soil is already saturated could trigger a mass 

18 movement. Studies have established relationships between antecedent precipitation and rainfall 

19 intensity thresholds for several different geographic regions (Glade et al. 2000; Godt et al. 2006; 

20 Chleborad et al. 2006). Moving forward, this research will consider the joint relationship 

21 between antecedent precipitation and precipitation intensity to better characterize potentially 

22 susceptible regions based on weekly, monthly or seasonal precipitation accumulation. 

23 
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1 Fig. 7 attempts to summarize the rainfall-landslide relationships over the three test areas. The top 

2 panel indicates that monthly rainfall is fairly well correlated with the exceedance index based on 

3 daily rainfall for the three test areas. Although one would think that extreme daily rainfall would 

4 be more closely associated with landslides, it is clear that the two rainfall statistics are related. 

5 The middle and bottom panels show that both the monthly rainfall and exceedance index are 

6 correlated to landslides, but that there is significant noise. In very approximate terms, a doubling 

7 of monthly rainfall from 150 to 300 mm is related to a fatal landslide increase of about a factor 

8 of three. A similar or slightly larger increase in fatal landslides is associated with a doubling of 

9 the exceedance index. These results are only indicative of areas that are already prone to 

10 landslides. Despite the data challenges intrinsic to this empirical approach, results shown here 

11 suggest that the GLC is very useful in estimating rainfall-landslide relations both in particular 

12 regions and even integrated over the globe. The data products evaluated here represent a very 

13 noisy process, particularly when evaluated across the globe. Despite this fact, we anticipate that 

14 if this evaluation were expanded to other study areas with sufficient numbers of landslide events, 

15 we may observe a more robust relationship between landslide reports and precipitation signals. 

16 5. Conclusions 

17 One of the unique aspects of the GLC is that it provides the first openly available, global picture 

18 of rainfall-triggered landslides over multiple years that can be compared with global precipitation 

19 estimates. Through the use of this catalog, the distribution and frequency of landslides and fatal 

20 landslides can be compared to distributions of satellite rainfall to better quantify these 

21 relationships. This analysis also allows us to evaluate the co-occurrence of extreme precipitation 

22 and landslide “hotspots” at large spatial scales and determine how landslide variations are related 

23 to meteorological changes. From analysis of the 2010 precipitation signatures over the three 
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1 study areas, it is clear that an observable signal exists between increases in reported and fatal 

2 landslide activity and increases in precipitation accumulation and daily intensity. It is not clear 

3 from the analyses and associated statistics that daily rainfall exceedance values are a better 

4 indicator of increased landslides than just simply the anomaly in monthly rainfall. The relative 

5 importance of daily extremes vs. monthly anomalies should be examined more thoroughly with 

6 additional data as the landslide catalog increases in length. Future analyses should also take into 

7 account regional or local differences in surface characteristics, such as are contained in landslide 

8 susceptibility indexes (e.g. Nadim et al. 2006; Guzzetti et al. 2005; Lepore et al. 2011). While 

9 other factors can modify this relationship including anthropogenic modification and tectonic 

10 weakening of hillslopes among others, understanding the relative distribution of extreme 

11 precipitation may help to shed new light on potential landslide activity at daily, monthly and 

12 yearly scales. We plan to re-evaluate these changes once we have built a larger record of 

13 reported events. 

14 

15 Through the type of study shown here, we may be able to better characterize the relative 

16 relationship between precipitation activity and potential landslide triggering and identify where 

17 landslides may impact populations based on natural variability in seasonal precipitation from 

18 teleconnections such as ENSO. Projections of precipitation intensity and distribution in a warmer 

19 world suggest that despite model uncertainties, rainfall in many of the monsoonal regions and 

20 tropical cyclone areas will likely become more extreme (IPCC 2007). One future direction of this 

21 study is to establish more concrete global relationships between extreme precipitation and 

22 landslide activity in order to better understand how landslide disasters may be modulated under 

23 climate change conditions. New satellite missions such as the Global Precipitation Measurement 
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1 (GPM) mission (www.gpm.nasa.gov) will also help to improve spatiotemporal coverage of 

2 precipitation measurements, enabling an extended record of satellite precipitation in order to 

3 better characterize the seasonal, yearly and decadal variability of extreme precipitation and its 

4 impact on landslide activity at the global scale. 

5 
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1 Figure Captions 

2 

3 Fig. 1 : Distribution of GLC for 2007 - 2010, highlighting: a) the distribution of total fatal 

4 landslides by country over the record with individual event locations for all reported landslides 

5 shown for 2010; and monthly distribution of b) fatal landslides and c) all reported landslides 

6 shown by year. 

7 

8 Fig. 2: Distribution of landslide reports for the years 2007 - 2010, showing a) reported and b) 

9 fatal landslides in North and South America, and c) reported and d) fatal landslides in Asia and 

10 Oceania. The boxes denote the three study areas evaluated in this paper: Central America, 

11 Himalayan Arc, and central-eastern China. Circles denote landslides for 2010, + signs display 

12 other years. The color denotes their moth of occurrence. 

13 

14 Fig. 3 : 2010 daily precipitation anomalies computed from a TMPA daily climatology for 1998 - 

15 2010 in mm/day. Blue (positive) areas indicate regions with higher daily precipitation totals, 

16 orange (negative) areas display dryer conditions for 2010. 

17 

18 Fig. 4: Precipitation analysis results for Central America study area, showing: a) Monthly 

19 rainfall accumulation for 2010 (red) with 12-year monthly climatology (green) calculated from 

20 the TPMA record (1998-2009); b) normalized threshold exceedance values (using the regional 

21 39 mm/day threshold) summed for each month in 2010 (red) and average values for 2007 - 2009 

22 (blue) compared to the landslide occurrence for 2010 and average number of reports from 2007 - 

23 2009; c) Q-Q plot showing the distribution of quantiles for the 12-year TMPA record (x-axis) vs. 

24 the 2010 daily values (y-axis). The interquartile line (red) and 1:1 line (green) provide a 
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1 reference to compare the distributions of quantiles for both periods. Evaluation statistics are 

2 shown in Table 1. 

3 

4 Fig. 5: Precipitation analysis results for Himalayan study area, showing: a) monthly climatology 

5 comparing 2010 (red) with 12-year climatology (green); b) normalized threshold exceedance 

6 values using the globally 79 mm/day threshold for 2010 and 2007 - 2009 with reported landslide 

7 events; c) Q-Q plot showing the distribution of quantiles for the 12-year TMPA record (x-axis) 

8 vs. the 2010 daily values (y-axis), compared against the 1:1 line (green) and interquartile line. 

9 

10 Fig. 6: Precipitation analysis results for the China study area, showing: a) monthly climatology 

11 comparing 2010 (red) with 12-year climatology (green); b) normalized threshold exceedance 

12 using the globally 79 mm/day threshold compared to landslides over the same periods; c) Q-Q 

13 plot showing the distribution of precipitation quantiles for the 12-year TMPA record (x-axis) vs. 

14 the 2010 daily values (y-axis). 

15 

16 Fig. 7: Scatter plots showing the monthly values for 2007 - 2010, comparing monthly 

17 precipitation and summed threshold exceedance pixels over each of the study areas with fatal 

18 landslides for each corresponding month over the 4 year record. The Central American region 

19 uses the 39 mm/day regional threshold (Guzzetti et al. 2008), while the Himalaya and China 

20 regions both use the 79 mm/day global threshold. Filled in symbols denote 2010 months. The 

21 three plots show: a) monthly rainfall (x-axis) vs. the sum of the exceedance values (y-axis); b) 

22 monthly rainfall (x-axis) vs. fatal landslides for each month (y-axis), showing the mean number 

23 of fatal landslides for each 50 mm; c) sum of exceedance values over each area (x-axis) vs. fatal 
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landslides (y-axis), showing the mean number of fatal landslides plotted at an interval of 50 


“hits” of the exceedance threshold. 


Table 1: Test statistics for the 3 study areas, showing the 75 th quantile, K-S test statistic and p- 


value, if the null was rejected, and the confidence level for rejecting the null. 






Number of Fatal Landslides 
(2007 - 2010 ) 

No landslide reports 
Reports but no fatal events 


] 6-10 
9 11-25 
| 26 - 50 
| 50+ 

Reported landslides in 2010 




Fig- 1 : Distribution of GLC for 2007 - 2010, highlighting: a) the distribution of total fatal 


landslides by country over the record with individual event locations for all reported landslides 


shown for 2010; and monthly distribution of b) fatal landslides and c) all reported landslides 


shown by year. 
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Fig. 2: Distribution of landslide reports for the years 2007 - 2010, showing a) reported and b) 


fatal landslides in North and South America, and c) reported and d) fatal landslides in Asia and 
Oceania. The boxes denote the three study areas evaluated in this paper: Central America, 
Himalayan Arc, and central-eastern China. Circles denote landslides for 2010, + signs display 
other years. The color denotes their moth of occurrence. 
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2 Fig. 3 : 2010 daily precipitation anomalies computed from a TMPA daily climatology for 1998 - 

3 2010 in mm/day. Blue (positive) areas indicate regions with higher daily precipitation totals. 
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2010 daily precip quantiles (mm) Number of reported landslides 


orange (negative) areas display dryer conditions for 2010. 
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Fig. 4: Precipitation analysis results for Central America study area, showing: a) Monthly 


rainfall accumulation for 2010 (red) with 12-year monthly climatology (green) calculated from 


1 the TPMA record (1998-2009); b) normalized threshold exceedance values (using the regional 

2 39 mm/day threshold) summed for each month in 2010 (red) and average values for 2007 - 2009 

3 (blue) compared to the landslide occurrence for 2010 and average number of reports from 2007 - 

4 2009; c) Q-Q plot showing the distribution of quantiles for the 1 2- year TMPA record (x-axis) vs. 

5 the 2010 daily values (y-axis). The interquartile line (red) and 1:1 line (green) provide a 

6 reference to compare the distributions of quantiles for both periods. Evaluation statistics are 

7 shown in Table 1. 
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2 Fig. 5: Precipitation analysis results for Himalayan study area, showing: a) monthly climatology 

3 comparing 2010 (red) with 12-year climatology (green); b) normalized threshold exceedance 
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1 values using the globally 79 mm/day threshold for 2010 and 2007 - 2009 with reported landslide 

2 events; c) Q-Q plot showing the distribution of quantiles for the 12-year TMPA record (x-axis) 

3 vs. the 2010 daily values (y-axis), compared against the 1:1 line (green) and interquartile line. 
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2 Fig. 6: Precipitation analysis results for the China study area, showing: a) monthly climatology 

3 comparing 2010 (red) with 12-year climatology (green); b) normali z ed threshold exceedance 

4 using the globally 79 mm/day threshold compared to landslides over the same periods; c) Q-Q 
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1 plot showing the distribution of precipitation quantiles for the 12-year TMPA record (x-axis) vs. 


2 the 2010 daily values (y-axis). 
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1 Fig. 7: Scatter plots showing the monthly values for 2007 - 2010, comparing monthly 

2 precipitation and summed threshold exceedance pixels over each of the study areas with fatal 

3 landslides for each corresponding month over the 4 year record. The Central American region 

4 uses the 39 mm/day regional threshold (Guzzetti et al. 2008), while the Himalaya and China 

5 regions both use the 79 mm/day global threshold. Filled in symbols denote 2010 months. The 

6 three plots show: a) monthly rainfall (x-axis) vs. the sum of the exceedance values (y-axis); b) 

7 monthly rainfall (x-axis) vs. fatal landslides for each month (y-axis), showing the mean number 

8 of fatal landslides for each 50 mm; c) sum of exceedance values over each area (x-axis) vs. fatal 

9 landslides (y-axis), showing the mean number of fatal landslides plotted at an interval of 50 

10 “hits” of the exceedance threshold. 
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75 th 

quantile 

(mm/day) 

K-S test stat 

p-value 

Null 

rejected 

alpha 

Central America 

6.63 

.1792 

0.0026 

Yes 

.997 

Himalayan Arc 

3.02 

.2119 

0.0004 

Yes 

.999 

China 

2.42 

0.0985 

0.3776 

No 

n/a 

China - 90 th 
Quantile 

5.97 

0.2213 

0.0553 

Yes 

.96 


12 

13 Table 1: Test statistics for the 3 study areas, showing the 75 th quantile, K-S test statistic and p- 

14 value, if the null was rejected, and the confidence level for rejecting the null. 
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